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Abstract 

We study steady streaming in a channel between two parallel permeable walls induced by 
oscillating (in time) blowing/suction at the walls. We obtain an asymptotic expansion of the 
solution of the Navier-Stokes equations in the limit when the amplitude of the normal displacements 
of fluid particles near the walls is much smaller that both the width of the channel and the thickness 
of the Stokes layer. It is demonstrated that the magnitude of the steady streaming is much bigger 
than the corresponding quantity in the case of the steady streaming produced by vibrations of 
impermeable boundaries. 

1 Introduction 

It is well-known that an oscillating (in time) body force or vibrations of the boundary of a domain 
occupied by a viscous fluid can produce not only an oscillating flow but al s o a (relative ly) weak steady 
flow, which is usually called steady streaming (see Lighthill . 1978 : Riley, 1967 . 200ll ). In this paper, 
we present a theory of steady streaming in a channel with fixed but permeable walls produced by 
given velocity at the walls which is oscillating in time with angular frequency uj. The basic parameters 
of the problem are the inverse Strouhal number and the inverse Reynolds number u, defined by 
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(1.1) 



where V^* is the amplitude of the oscillating velocity at the walls, d is the distance between the walls 
and u* is the kinematic viscosity of the fluid. Parameter measures the ratio of the amplitude of 
the displacements of fluid particles in an oscillating velocity field with amplitude Vq to the distance 
between the walls. Another dimensionless parameter which is widely used in literature is the 'streaming 



Reynolds number' Rg 



In terms of parameters e and u, Rg = e jv. We are interested in 



the asymptotic behaviour of solutions of the Navier-Stokes equations in the limit e ^ 1 and v ^ \. 
This means that the amplitude of displacements of fluid particles is much smaller than the thickness 
of the Stokes layer. Note that this limit corresponds to small Rg {Rg ~ <C 1). 

Early studies of the steady streaming in a channel induced by vibrations of the walls had been 
focused on the problem of peristaltic pumping in channels and pipes under the assumption of low 
Reynolds numbers {R < C 1) and small amplitude-to- wavelength ratio (see, e.g., IJafFrin &: Shapiro! . 



197ll : IWilson Pantonl . Il979l ). In recent years, there had been considerable rene wed interest in the 



prob l em motivated by poss i ble applications of steady strea ming to micro-mixing (jSelverov &: Stond . 



200ll: lYi. Bau. h Hul. 20021: ICarlsson. Sen h Lofdahll . |2005| ) and to drag reduction in channel flows 



dHoepffner fc Fukagatal . l2009l ~ In all asymptotic theories of the steady streaming produced by vi- 
brating impermeable boundaries, the magnitude of steady streaming is O(e^) for e <C 1. This is 
because steady streaming results from the nonlinearity in the Navier-Stokes equations and parameter 
measures the magnitude of the nonlinear term. The aim of the present study is to show that if the 
boundary is permeable, then the magnitude of steady streaming is 0(e) for small e, i.e. much bigger 
than in the case of an impermeable boundary. 

To obtain an as ymptotic expansion, we emp loy the Vishik-Lyus ternik method (see, e.g. 
197nl : iNavfehl . Il973l ). It is described in detail in lllin fc Sadkj (I2OI0I I. 



Trenoginl . 
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The outhne of the paper is as fohows. In Section 2, we formulate the mathematical problem. 
In Section 3, the asymptotic expansion of the solution is described. In Section 4, we consider sim- 
ple examples in which the leading order asymptotic solution can be obtained analytically. Finally, 
conclusions are presented in Section 5. 



2 Formulation of the problem 

We consider a two-dimensional viscous incompressible flow in an infinite channel of width d. The walls 
of the channel are permeable for the fluid, and the flow is produced by a given velocity at the walls 
which is assumed to be periodic along the channel with period L* and oscillating in time with angular 
frequency uj. We will use the following non-dimensional quantities: 

X* V* p* 

T = u;t, x=— , v = — , p- 



d ^0 P^^^o* 

Here t* is time; x* = (x*,y*); x* and y* are Cartesian coordinates, the x* axis being parallel to 
the channel; v* = (n*,v*) is the velocity of the fluid; p* is the pressure; p is the density; Vq* is the 
maximum of the given velocity at the walls over all x* and t* . In these variables, the Navier-Stokes 
equations take the form 

v^ + e^(v- V)v = -Vp + e^i/VV, V • v = 0, (2.1) 

where the dimensionless parameters and u are defined by Eq. (jl.ip . Equations (|2.ip are to be 
solved subject to the boundary conditions 

v|j,=o = ^^(x, r, e), w\y=i = Y\x, r, e). (2.2) 

Here = {U"", V"") and = (C/^, V^) are given functions which are 27r-periodic in r and have zero 
mean value: 

V" = — /"v"(x,T)dr = 0, Y^ = — [Y\x,T)dT = ^. (2.3) 
271 J 2-K J 



They are also periodic in x with period L = L* /d and satisfy the condition 

L L 

V''{x,T)dx = [ V\x,T)dx (2.4) 



which follows from incompressibility of the fluid. In what follows we are interested in the asymptotic 
behaviour of periodic (both in r and x) solutions of Eqs. (j2.ip . (j2.2p in the limit e — )• and = 0(1). 
We assume that V'*'^(x,r, e) can be presented in the form 

V"'^(x, T, 6) = Vf{x, t) + 6 (x, r) + e^vf (x, r) + . . . (2.5) 

We seek a solution of (12. ip , (j2.2p in the form 

u = u'' + u'' + vl', V = v"" + ev" + ev^, p = p"" + p" + p^ . (2.6) 

Here vJ" , , p^ are functions of x, y, r and e; n", depend on x, r, e and the boundary layer 

variable ^ = y/e; u'', v'', depend on x, r, e and the boundary layer variable rj = {1 — y) / e. Functions 
u^, v^, p^ represent a regular expansion of the solution in power series in e (an outer solution), and 
(li*^, v", p") and (u'', , p^) correspond to boundary layer corrections to this regular expansion (inner 
solutions). We assume that the boundary layer parts of the expansion rapidly decay outside thin 
boundary layers: 

u^'^v^-^p"- = o{i~'') as i^oo and ,p^ = oir]'") as r] ^ oo (2.7) 

for every s > 0. This assumption will be verified a posteriori. 
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3 Asymptotic expansion 

3.1 Regular part of the expansion 

Let 

V =V^ + ewl + e\l + =pl + ep[ + e^p'2 + (3.1) 

where = {vJ',v'^), vjl = (n^jV^) {k = 0,1,2,...). The successive approximations v^, p^^ {k = 
0, 1, 2, . . . ) satisfy the equations: 

drV^ = -Vpl + Yk, V-v^, = 0, (3.2) 

where Fq = 0, Fi = and 

fe-2 

= - E(< • ^)^l-2-i + ^VVL2 (3.3) 

1=0 

for k > 2. In what fohows, we will use the following notation: for any 27r-periodic /(r), 

f{r)=f + f(.r), f= ^ I f{r)dT (3.4) 



where / is the mean value of /(r) and, by definition, /(r) = /(r) — /. 
3.1.1 Leading-order equations 

A general solution of Eqs. ()3.2p for /c = 0, which is periodic in r, can be written as 

v5 = vg + vo, Vq = "^4)0 (3-5) 
where <pQ has zero mean value and is the solution of the boundary value problem 

V^(^o = 0, (l}oy\y=o = VQ{x,T), (l)Qy\y=i = V^{x,T), (t)o{x + L , y) = (t)Q{x , y) . (3.6) 



The boundary conditions at y = and y = 1 in (j3.6p will be justified later. 

On averaging the equation for ^^2 (the first equation (13. 2p for k = 2) and the incompressibility 
condition for Vg and using the fact that Vq is irrotational, we obtain 

(v5-v)v5 = -vno + i^v2v5, v-v5 = o, (3.7) 



where XIq = P2 + l^'/'oP/2. Equations (|3.7p represent the time-independent Navier-Stokes equations. 
It will be shown later that boundary conditions for Vq are 

^l\y=Q = ^l\y=l = 0. (3.8) 

The only solution of (|3.7p that is periodic in x and satisfies p.Sp is zero solution: 

v5 = 0. (3.9) 
This means that there is no steady streaming in the leading order of the expansion. 
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3.1.2 First-order equations 

The solution of Eqs. (j3.2p for k = 1 can be written as 

v^ = v^ + vi, v^ = V0i (3.10) 
where 0i has zero mean value and is the solution of the boundary value problem 

V^(/)i = 0, (piy\y=o = ai{x,T), (piy\y=i = bi{x,T), (j)i{x + L , y) = (j)i{x , y) . (3.11) 

Functions ai(a;,r) and 6i(x,r) will be defined later. 

To obtain equations for v'[{x.), we average the incompressibility condition for and the equation 
for Vg (the first equation ()3.2p for k = 3) and then take account of Eq. (13.90 and the fact that Vq and 

are both irrotational. This yields 

= -VHi + z/V^v^, V • = 0. (3.12) 

where Hi = + {Vcpo • Vcpi). Thus, the first order averaged outer flow is described by the Stokes 
equations. Boundary conditions for will be specified later. 

3.2 Boundary layers near the walls at y = and y = 1 
3.2.1 Boundary layer equations 

To derive boundary layer equations near the bottom wall, we ignore u'^, and p^, because they are 
supposed to be small everywhere except a thin boundary layer near y = 1, and assume that 

u = u''Q + ul + e{ul+u'{) + v = vl + e{v\ + vl) + p = p'^ + pi + e{p\ + pi) + . . . (3.13) 

We substitute these into Eq. (|2.ip and take into account that u^, f^, satisfy (j3.2p . Then we make 
the change of variables y = e^-, expand every function of in Taylor's series at e = and collect 
terms of the equal powers in e. As a result, we obtain 

<r+Pl.-^<ii = Fk^ Pk = Gt, ut, + vt^ = (3.14) 

for A; = 0, 1, . . . In Eqs. p.l4p . functions and depend on Vq, . . . , v^_j^, Uq, . . . , u^^^, Vq, . . . , v^^^. 
For k = 0, 1, these are given by 

FS = 0, Ff = -yo"(x,r)<^, Gg = 0, = 0. (3.15) 

A similar procedure leads to the equations of the boundary layer near the upper wall: 

''At + Pkx ~ '^''Arjr] = Pkr] = ""fca; ~ ''^kr] = (3.16) 

for A; = 0, 1, ... . Functions F^ and for A; = 0, 1 are given by 

Fo^ = 0, F^ = Vo\x,t)4,,, Gl = 0, G5 = 0. (3.17) 

Before we describe how the above boundary layer equations can be solved, we need to discuss boundary 
conditions that must be satisfied by terms of each order in our asymptotic expansion. 
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3.2.2 Boundary conditions 

In view of (j2.7p . we require tliat for every s > and for each = 0, 1, . . . , 



'^k^vl,pl = o{i as i^oo and 4, Pfc = o(r/ ^) as r/ ^ oo. 



(3.18) 



Now we substitute (|3.13p into the first Eq. (12. 2p and collect terms of equal powers in e. This leads to 
the following boundary conditions at y = 0: 



Ui. 



y=0 



y=0 



5=0 



USix,T), Vl 



y=0 



+ v 



5=0 



for A; > 1. Similarly, boundary conditions at y = 1 are given by 



y=l 



+ 4 



y=l 



+4 



r?=0 



7?=0 



y=l 



y=l 



+4-1 



77=0 



for k > 1. Note that (|3.19p and (|3.2ip justify boundary conditions for (j)Qy in (|3.6 



(3.19) 
(3.20) 

(3.21) 
(3.22) 



3.2.3 Leading order equations 



Boundary layer at y = 0. In the leading order, the boundary layer at y = is described by (|3.14p for 
k = 0. The condition of decay at infinity (in variable ^) for and the second Eq. (j3.14p imply that 
Pq = 0. Hence, the first Eq. (I3.14p reduces to the heat equation 



Ot 



Boundary condition for at ^ = follows from ()3.19p and is given by 



u, 



oL=o= US{x,t) -(l)Oa,{x,0,T). 



(3.23) 



(3.24) 



The periodic (in r) solution of Eq. (I3.23P must satisfy (I3.24p and the condition of decay at infinity. 

On averaging Eq. (|3.23p . we find that ng^^ = 0. The only solution of this equation that satisfies 
the decay condition at infinity and the boundary condition ^01^=0 = (which follows from (I3.24p and 
the fact that Uq = and (j)Q = 0) is zero solution. Thus, in the leading order the boundary layer at 
y = is purely oscillatory: = 0. 

The normal velocity Vq is determined from the third Eq. (j3.14p for k = 0: 



V^ix,(,T)=d, / txg(x,C',T)dC'. 



(3.25) 



Here the constant of integration is chosen so as to guarantee that — )> as ^ — )• oo. According 
(|3.13p . Vq{x,^,t) enters the first-order term in the expansion and, therefore, Vq gives us the 
boundary condition for the next approximation of the outer solution. Indeed, in view of (j3.20p (with 
k = 1), we have v\{x,y,T) \y^Q= ^i{x,t) — Vq{x,0,t), so that function ai(x,T) in problem (I3.1ip is 



ai{x,T) = V^{x,t) 



Vq(X, 



0,t). Note also that Eq. ([3:25|) imphes that v^ = 0. 



Boundary layer at y = 1. Exactly the same arguments as above lead to the problem 



4r 
4' 



Tb 



r7=0 



Uq{x,t) - (t)Qx{x,l,T), tig as 7?-^oo. 



(3.26) 
(3.27) 
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Again, it follows from (j3.26p and (|3.27p that Uq = 0. The normal velocity is given by 

oo 

v^oix, r], t) = -d^ [ 4ix, rj', T)dr]' . (3.28) 



Then it follows from (13:22]) (with A: = 1) that h in problem (l3TT|) is bi{x,T) = Vf{x,T) - v^{x,0,t). 
Again, K28\i implies that = 0. 

Averaged outer flow. On averaging (j3.19p and (|3.2ip and using the fact that both boundary layers 
are purely oscillatory and our assumption that Vq'^ = 0, we arrive at conclusion that Vq |y_o= 

and Vq |y_^= 0- This justifies boundary conditions (13. 8p and our conclusion that there is no steady 
streaming in the leading order of the expansion. 

3.2.4 First-order equations 

Oscillatory outer flow. Since functions ai{x,T) and bi{x,T), which appear in problem (|3.1ip . are now 
known, the problem can be solved, thus giving us v^(x,y,r). 

Boundary layer at y = 0. Consider now Eqs. (j3.14p for k = 1. Again, the condition of decay at 
infinity for pi and the second Eq. (|3.14p imply that pi = 0. Hence, we have 

< = z.n?g^-yj^(x,r)7/S^. (3.29) 



Averaging this equation, we find that I'u^^ = Vq{x,t)uq^. Integration over ^ yields 



V^-{x,T)ul{x,i',T)di'. (3.30) 



Here the constants of integration are chosen so as to satisfy the condition of decay at infinity. According 
to p.20p . we have 



oo 

1 



ul\y=o = -ul\^=o = -J Vo^{x,T)u^{x,C',T)dC', (3.31) 



The oscillatory part of uf, as well as both averaged and oscillatory parts of vf can also be found but 
are not needed in what follows. 

Boundary layer at y = 1. A similar analysis leads to 

oo 

1 



I' 



Vo'{x,T)u'oix,r,',T)dr,', (3.33) 



oo 

1 



h\y=i- - 



Vo'{x,T)uUx,v',r)dr,', (3.34) 



vl\y^i=0- (3.35) 

Averaged outer flow. The first order averaged outer flow is described by the Stokes equations (j3.12p 
subject to the boundary conditions (I3.3ip . (13.321) . (13.341) and (I3.35p . This boundary value problem 
describes a steady Stokes flow produced by a given distribution of the velocity at the boundary. It will 
have a non-zero solution provided that at least one of the integrals on the right sides of Eqs. (13.3ip 
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and (j3.34p is non-zero. In this case, steady streaming is the effect of first order in e. This is in contrast 
with steady streaming produced by vibrating impermeable wahs where a non-zero averaged velocity 
appears in the second order of the expansion. Thus, the steady component of the flow in our problem 
is much stronger than that in the case of vibrating impermeable walls. 

Leading order asymptotic for the averaged flow. To summarise, we have obtained the first non-zero 
term in the expansion of the steady component of the flow. The averaged velocity field has the form 

u = e{ul + u1 + u\) +0{e^), v = evl + 0{e^) 

where boundary layer contributions uf and u\ are given by Eqs. (j3.30p and (j3.33p and where is the 
solution of the Stokes equations that satisfies boundary conditions (j3.3ip . (|3.32p . (|3.34p and (|3.35p . 
If we introduce the stream function for the averaged flow ^jJ defined by the standard relations u = ipy 
and V = —ipx, then the corresponding expansion of ip will have the form 

where ipl is the stream function for = (u^, -(;[). Note that boundary layer terms do not appear in 
the leading order of the expansion for ^. 

We note also that, in our problem, there is no need to consider the Stokes drift as it is the effect 
of higher order in e. 



4 Examples 

Below we consider a few simple examples in which the velocity at the walls oscillates harmonically 
in time. In all the examples below, the tangent velocity at the walls is zero = = 0), and 
the normal velocity is non-zero {V"" 7^ and 7^ 0). Examples with zero normal velocity and 
non-zero tangent velocity are not considered, as in this case, steady streami ng appears in highe r order 
approximations (in the case of half space this problem had been treated bv Iviadimirovl (ji^)). 



4.1 Example 1: Standing waves 

Let = cosA:3; cosre^ and V* = aV" where k = 2-k/L and a = ±1. This choice corresponds to 
standing waves of injection/suction applied at the boundaries of the channel (a = 1 if the waves are 
in phase, and a = — 1 if they have opposite phase). After substitution of V" and V'' in the general 
formulae of Section 3, we find that 

-ri -ri ^(^) ■ /r,, ^ An\ COsh{k) - a 

A\/2v smh(/c) 
Solving the Stokes equations (j3.12p with boundary conditions ()3.32p . ()3.35p and (|4.ip . we obtain 



4- 



1 A(k) 

7^ sinh(2fc) - 2k sinh(2%) + y sinh [2fc(l - y)]} sin(2fcx). (4.2) 

Typical streamlines are shown in Fig. (U^a). It is clear from ()4.2p that the flow picture is the same 
for a = 1 and a = — 1, the only difference is in the magnitude of the flow. The latter is determined 
by Aifi). For a = 1, ^ is an increasing function, A{k) ~ k for k \ and Aifi) — )• 1 as A; — )• 00. For 
a = —1, A is a decreasing function, A{k) ~ 1/k for A; ^ 1 and A(k) — )• 1 as A; — )• 00. So, for all 
k > 0, the magnitude of the steady streaming for a = — 1 is greater than for a = 1. Thus, the steady 
streaming is stronger when standing waves of suction/injection applied at the walls have opposite 
phase. Also, in this case the magnitude of the steady streaming increases with the wavelength of the 
waves. 
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Figure 1: The streamlines ipl = const for u = 1, L = 3 and a = 1: (a) standing waves; (b) waves 
travelling in opposite directions. 



4.2 Example 2: Waves travelling in the same direction 

Let now the velocity at both walls be purely normal and have the form of waves travelling in the 
direction of the x axis: V" = cos{kx — T)ey and = aV" where k and a are the same as in 
Example 1. Boundary conditions (j3.3ip and (j3.34p take the form 

A{k) 



2V2z^ 



(4.3) 



where A{k) is given by ()4.ip . The Stokes equations (j3.12p subject to (j3.32p . (j3.35p and (j4.3p lead to the 
constant solutiorH: u\ = —A{k)/2^/2i/, vl = 0. Thus, the waves travelling in the same direction pro- 
duce a constant mean flow whose direction is opposite to the direction in which t he waves advance. This 



(somewhat surprising) result agrees with numerical simulations reported in (jHoepffner &: Fukagatal . 



20091 ). The magnitude of the mean flow is determined by A(k). Properties of A{k) imply that the 



most efficient way to generate the mean unidirectional flow is to apply suction/blowing in the form of 
waves travelling in the same direction and having opposite phase. 



4.3 Example 3: Waves travelling in the opposite directions 



Let now the normal velocity at the walls have the form of waves travelling in opposite directions: 
V" = cos(A:x — r) e^^ and = a cos{kx + r) where k and a are the same as in Examples 1 and 2. 
Boundary conditions for reduce to 



*lly=0 



1 B-{k,x) 



2\/2i^ sinh(A;) 



Ul\y=l 



1 B+{k,x) 



2V2u sinh(A;) 



(4.4) 



where = cosh(A;) + a[cos(2A;x) it sin(2A;x)]. The corresponding solution of the Stokes equations is 
given by 



1 



2V2i' 



cosh(/c) 
sinh(A;) 



y(l -y) + D-y sinh[2/c(l - y)] + D+{1- y) sinh(2%) 



(4.5) 



^The Stokes equations (|3.12p also admit solutions with a nonzero pressure gradient VIIi = cqQx (cq — const), which 
are not considered here, because this would be equivalent to a modification of our problem, allowing the presence of a 
weak O(e^) pressure gradient. 
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smn(A;j 



where 

cos(2A;x) sm.{2kx) 
sinh(2/c) + 2k sinh(2/c) - 2k _ 

Typical streamlines of the flow (|4.5p for a = 1 are shown in Fig. [Hb). In the case of a = —1, the flow 
picture is the same except that it is shifted along x axis by L/4 (this can be deduced directly from 

CSD). 



5 Conclusions 

We have considered incompressible flows in a channel between two parallel permeable walls and con- 
structed an asymptotic expansion of solutions of the Navier-Stokes equations in the limit when the 
amplitude of displacements of fluid particles near the walls is mich smaller than both the width of 
the channel and the thickness of the Stokes layer. The asymptotic procedure is based on the Vishik- 
Lyusternik method and can be used to construct as many terms of the expansion as necessary. In 
the leading order, the averaged flow is described by the stationary Stokes equations subject to the 
boundary conditions that are determined by the boundary layers at the walls. The key difference 
between the present expansion and the asymptotic theories of steady streaming induced by vibrating 
impermeable boundaries is that, in our study, the magnitude of the steady streaming is 0(e), which 
is much bigger than O(e^) steady streaming in the case of impermeable walls. 

The general formulae have been applied to three particular examples of steady streaming induced 
by blowing/suction at the walls in the form of standing and travelling plane waves. In the case of 
standing waves, the averaged flow has the form of a double array of vortices (see Fig. HJa)). For 
short waves the vortices are concentrated near the walls, while long waves produce vortices that fill 
the entire channel. 

If the normal velocity at the walls have the form of plane harmonic waves which travel in the 
same direction, the induced steady flow is a constant unidirectional flow whose direction is opposite 
to the direction in which the waves travel. This is different from the case of the steady streaming 
generated by vibrations of impermeable walls where the induced flow had the same direction as the 
travelling wave. As far as we are aware, this was first observed in numerical simulations performed by 
Hoepffner fc Fukagatal jioO^) . 



If the normal velocities at the walls have the form of plane waves traveling in opposite directions, 
the averaged flow is a superposition of a shear flow with a linear velocity profile and a periodic array 
of vortices ('cat's eyes') in the center of the channel. When the wavelength is small, the vortices are 
weak. Their intensity and size monotonically grow with the wavelength and eventually they fill the 
entire gap between the walls. 

There are many open problems in this area. In particular, it is not clear how the present theory 
can be extended to the case of Rs ~ 1. Although the problem does not involve a moving boundary 
which, in general, simplifies things, there is a technical difficulty of a different sort. It is related to the 
leading order boundary layer equations for ~ 1 which are difficult to solve ana lyticalljU. This is a 
subject of a continuing investigation. 

Acknowledgments. I am grateful to A. B. Morgulis and V. A. Vladimirov for stimulating discussions. 
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